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The total r hadronic width can be accurately calculated using analyticity and 
the operator product expansion. The theoretical analysis of this observable is up¬ 
dated to include all available perturbative and non-perturbative corrections. The 
experimental determination of (y.s{M^) and its actual uncertainties are discussed. 


1 Introduction 

The inclusive character of thej-total r hadronic width renders possible an ac- 
curate calculation of the ratioEra [( 7 ) represents additional photons or lepton 
pairs] 

_ r[T- ^ hadrons ( 7 )] 

At — / \1 ) O 

r[r ^ Vre VeVi)\ 

using standard field theoretic methods. If strong and electroweak radiative 
corrections are ignored and if the masses of final-state particles are neglected, 
the universality of the W coupling to the fermionic charged currents implies 

i?, ~ 3, (2) 

which compares quite well with the experimental average0 Rt = 3.649±0.014. 
This provides strong evidence for the colour degree of freedom Nc- 

The QCD dynamics is able to account quantitatively for the 20% difference 
between the naive prediction (H) and the measured value of Rr- Moreover, the 
uncertainties in the theoretical calculation of Rt are quite small. The value of 
Rt can then be accurately predicted as a function of as{M^). Alternatively, 
measurements of inclusive r decay rates can be used to determine the value 
of the QCD running coupling as{M^) at the scale of the r mass. In fact, r 
decay is probably the lowest-energy process from which the running coupling 
constant can be extracted cleanly, without hopeless complications from non- 
perturbative effects. The r massjij = (1777.00 ilio' 27 ) MeV, lies fortuitously 
in a compromise region where the coupling constant as(M^) is large enough 
that Rt is sensitive to its value, yet still small enough that the perturbative 
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expansion still converges well. Moreover, the non-perturbative contributions 
to the total r hadronic width are very small. 

It is the inclusive nature of the total semihadronic decay rate that makes a 
rigorous theoretical calculation of Rr possible. The only separate contributions 
to Rr that can be calculated are those associated with specific quark currents. 
We can calculate the non-strange and strange contributions to Rr, and resolve 
these further into vector and axial-vector contributions. Since strange decays 
cannot be resolved experimentally into vector and axial-vector contributions, 
we will decompose our predictions for Rr into only three categories: 

Rr = Rr,v + Rt.A + Rt,S ■ ( 3 ) 

Non-strange semihadronic decays of the r are resolved experimentally into 
vector (Rr^v) and axial-vector contributions according to whether the 

hadronic final state includes an even or odd number of pions. Strange decays 
{Rr,s) are of course identified by the presence of an odd number of kaons in 
the final state. The naive predictions for these three ratios are i?r,y — Rt,a — 
(iVc/2)|I4dP and Rr^s - fVdKsP, which add up to (|). 


2 Theoretical Framework 


The theoretical analysis of Rr involves the two-point correlation functions for 
the vector and axial-vector colour-singlet 

quark currents {i,j = u, d, s): 


n^yiq) = 11 d^xe^‘^^{0\T{V:;{x)Vr^i0n0), (4) 

^ z I d^xe^<^^{0\TiA>^ix)A’^^iono}. (5) 

They have the Lorentz decompositions 

Kv/Ai<l) = ^^A^q^) + q^'^ ( 6 ) 

where the superscript (J) denotes the angular momentum J = 1 or J = 0 in 
the hadronic rest frame. 

The imaginary parts of the two-point functions are propor¬ 

tional to the spectral functions for hadrons with the corresponding quantum 
numbers. The semihadronic decay rate of the r can be written as an integral 
of these spectral functions over the invariant mass s of the final-state hadrons: 


Rr = 127r 



ds 




Imn*-^^ (s) -I- Imll^®^ (s) 


■ ( 7 ) 
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Figure 1: Integration contour in the complex s plane, used to obtain Eq. (g) 


The appropriate combinations of correlators are 

n(^)(s) ^ iKdl" (n2^(s) + n2^(s)) + (ngvW + • (8) 

The contributions coming from the first two terms correspond to Rt,v and 
Rt,a respectively, while i?T,s contains the remaining Cabibbo-suppressed con¬ 
tributions. 

Since the hadronic spectral functions are sensitive to the non-perturbative 
effects of QCD that bind quarks into hadrons, the integrand in Eq. ( 0 ) can¬ 
not be calculated at present from QCD. Nevertheless the integral itself can be 
calculated systematically by exploiting the analytic properties of the correla¬ 
tors They are analytic functions of s except along the positive real 

s axis, where their imaginary parts have discontinuities. The integral can 
therefore be expressed as a contour integral in the complex s plane running 
counter-clockwise around the circle |s| = M^: 


Rr = 67ri 


|=M2 


ds 

M2 






The advantage of expression (i) over (0) is that it requires the correlators 
only for complex s of order M^, which is significantly larger than the scale as¬ 
sociated with non-perturbative effects in QCD. The short-distance Operator 
Product Expansion (OPE) can therefore be used to organize the perturba- 
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tive and ntp-perturbative contributions to the correlators into a systematic 
expansiorJld in powers of 1/s, 


D=2n dimO=D ' 


( 10 ) 


where the inner sum is over local gauge-invariant scalar operators of dimension 
= 0, 2,4,... The possible uncertainties associated with the use of the OPE 
near the time-like axis are absent in this case, because the integrand in Eq. 
includes a factor (1 — s/M^)^, which provides a double zero at s = M^, effec¬ 
tively suppressing the contribution from the region near the branch cut. The 
parameter /i in Eq. ( 0 ) is an arbitrary factorization scale, which separates 
long-distance non-perturbative effects, which are absorbed into the vacuum 
matrix elements {0{^)), from short-distance effects, which belong in the Wil¬ 
son coefficients C^'^^(s,/i). The D = 0 term (unit operator) corresponds to the 
pure perturbative contributions, neglecting quark masses. The leading quark- 
mass corrections generate the D = 2 term. The first dynamical operators 
involving non-perturbative physics appear at I? = 4. Inserting the functions 
(^ into (||) and evaluating the contour integral, Rr can be expressed as an ex¬ 
pansion in powers of 1/M./, with coefficients that depend only logarithmically 
on Mr- 

It is convenient to express the corrections to Rr from dimension-H op¬ 
erators in terms of the fractional corrections , to the naive contribution 
from the current with quantum numbers zj, Y or zj, A: 


Rt,V/A 


Rt,S 



( 11 ) 


= {SifY+S^^\)/2 is the average of the vector and axial-vector corrections. 
The dimension-0 contribution is the purely perturbative correction neglecting 
quark masses, which is the same for all the components of Rr'. ~ 

The factorslii 



Sew 


fa^\ - f a{M^) \ - f a{Ml) \ - 
\a{M?)) \a{M^)) 


1.0194 



a(M/) 


0 . 0010 , 


( 12 ) 

(13) 
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contain the known electroweak corrections. 

Adding the three terms, the total ratio Rr is 

A. = 3 {\Vud? + \Vus?) Sew 1^1 + S'ew + 

+ +sin2 6»c<5if^)|, (14) 

where sin^ dc = |KsP/(|KdP + |KsP)- 


3 Perturbative Corrections 


In the chiral limit (m„ = rud = rus = 0), the vector and axial-vector cur¬ 
rents are conserved. This implies (fi) C therefore, only the correlator 
n(o-i-i) contributes to Eq. (^). Owingl^’llj to the chiral invariance of massless 
QCD, n|°+^^(s) = = n(s) (z ^ j) at any finite order in as- 

The result is more conveniently expressed in terms of the logarithmic 
derivative of the two-point correlation function of the vector (axial) current, 


D{s) 



1 

47r^ 


n—O 



(15) 


which satisfies an hpuMgeneous renormalization-group equation. The Kn co¬ 
efficients are knownt4E3 to order a^. For ny = 3 quark flavours, one has: 


Ko = if 1 = 1 : 


299 

if2 = ^-9C(3) = 1.63982, 

Kz{MS) = ^ - ^C(3) + yC(5) = 6.37101. 

The perturbative component of Rr is given by 

Rpert 


{\Vud? + \Vus?)Sew 
where the functionsi 


= 3{l + d(°)} = 3^if„A(”)(a,), 




ds f q:s(— s) 


|s|=M2 S 


n—0 


s s'* 


(16) 


(17) 


( 18 ) 
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are contour integrals in the complex plane, which only depend on as{M^). 

The running coupling as(—s) in Eq. ( p^ can be expanded in powers of 
as{M^), with coefficients that are polynomials in log {—s/M^). The pert^ba- 
tive expansion of in powers of ar = asiM^)/^ then takes the formiErlD 


= Y,{Kn+9n)<, 


n—1 

= CLt H" 


K2 - g/3i ) a? 


(19) 


19 19 


265 - 247r2 


f5l\al+0{at) 


288 

= ar + 5.2023 + 26.366 a? + O(a^). 

The complex integration along the circle |s| = generates the (?„ coeffi¬ 
cients, which depend on Km<n and on firn<n, where 


/3i = -9/2, /32 = -8, /33(M5) =-3863/192, ... (20) 

are the coefficients of the QCD (3 function for n/ = 3. One observesi that 
the gn contributions are larger than the original Kn coefficients (92 = 3.5625, 
93 = 19.9949). For instance, the bold-guess value K 4 it's(i^a 7 /^ 2 ) ~ 25 
is to be compared with 34 = 78.0029. Tl^p. large corrections give rise to a 
sizeable renormalization scale dependences’^ of the truncated O(a^) result. 
The reason of such uncomfortably large contributionsB stems from the long 
running along the circle s = exp (icp) (^e[0, 27r]) in Eq. (p^). When the 
running coupling ««(—s) is expanded in powers of as(M^), one gets imaginary 
logarithms, log (—s/M/) = i{(j) — tt), which are large in some parts of the 
integration range. The radius of convergence of this expansion is actually 
quite small. A numerical analysis of the seriesO shows that, at the three-loop 
level, an upper estimate for the convergence radius is aT,conv < 0.11. 

Note, however, that there is no deep reason to stop the integral 

expansions at 0{a^). One can calculate the expansion to all orders 

in as, apart from the unknown /3n>3 contributions, which are likely to be neg¬ 
ligible. Even for larger than the radius of convergence Or.conv, the integrals 
(as) are well-defined functions that can be numerically computed, by using 
in Eq. (|l8| ) the exact solution for as(—s) obtained from the renormalization- 
group /3-function equation. Thus a more appropriate approachQ’cd is to use a 
Kn expansion of as in Eq. (p?]), and to fully keep the known three-loop 
calculation of the functions A("^(as). The perturbative uncertainties are then 
reduced to the corrections coming from the unknown Pn >3 and Kn >3 contri¬ 
butions, since the gn contributions are properly resummed to all orders. To 
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appreciate the size of the effect, Table gives the exact resultfor A^'^\as) 
{n = 1,2,3) obtained at the one-, two- and three-loop approximations (i.e. 
/3„>i = 0, /3n>2 = 0, and /3„>3 = 0, respectively), together with the final value 
of for Qr =0.1. For comparison, the numbers coming from the truncated 
expressions at order are also given. Although the difference between the 
exact and truncated results represents a tiny 0.7% effect on Rr, it produces a 
sizeable 4% shift on the value of The (5^°^ shift, which reflects into a cor¬ 
responding shift in the experimental as{M^) determination, depends strongly 
on the value of the coupling constant; for Ur = 0.14 the shift reaches the 
—20% level. 


Table 1: Exact results for (n = 1,2,3) obtained at the one- (/3n>l = 0), two- 

iPn>2 = 0) and three—loop {j3n>z = 0) approximations, together with the final value of (5^^^, 
for ttr = 0.1. For comparison, the numbers coming from the truncated expressions at order 
are also given. 


Loops 





1 

0.13247 

0.01570 

0.00170 

0.1690 

2 

0.13523 

0.01575 

0.00163 

0.1714 

3 

0.13540 

0.01565 

0.00160 

0.1712 

“oRT 

0.14394 

0.01713 

0.00100 

0.1784 


Notice that the difference between using the one- or two-loop approxima- 
tion to the /3 function is already quite small (1.4% effect on (5^°^), while the 
change induced by the three-loop corrections is completely negligible (0.1%). 
Therefore (unless the (3 function has some unexpected pathological behaviour 
at higher orders), the error induced by the truncation of the (3 function at third 
order should be smaller than 0.1% and therefore can be safely neglected. The 
sensitivity on thexhaice of renormalization scale and renormalization scheme 
is also very smallBE3 

The dominant perturbative uncertainties come from the unknowiL-higher- 
order coefficients iF„> 3 . The 0{a'^) contribution has been estimatecupjsing 
scheme-invariant methods, namely the principle_pf minimal sensitivityEj and 
the effective charge approachwith the results 

ATf * = 27.5 . (21) 

This number is very close to the naive gV^S ~ {K^/K 2 )K 3 « 25. A 
similar estimate, = 24.8, is obtainedcZra in the limit of a large numbje:; 

of quark flavours, using the so-called naive non-abelianization prescriptions 
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Table 2: 5^^'^ for different values of as(M^). 




A((5W) 

0 

II 

Ki = 27.5 

0.30 

0.161 

0.164 

±0.006 

0.31 

0.168 

0.172 

±0.007 

0.32 

0.176 

0.180 

±0.008 

0.33 

0.183 

0.188 

±0.008 

0.34 

0.191 

0.196 

±0.009 

0.35 

0.198 

0.203 

±0.010 

0.36 

0.205 

0.211 

±0.010 

0.37 

0.213 

0.219 

±0.011 

0.38 

0.220 

0.226 

±0.012 

0.39 

0.227 

0.234 

±0.012 

0.40 

0.234 

0.241 

±0.013 


{rif ^ 3/3i = Uf — ^ = — Fron^ fit to the experimental t data, the value 
/Ff ‘ = 29 ± 5 has been also quotedS 

Using the estimate (^), the 0{a‘^) correction amounts to a 0.004 increase 
of for ar = 0.1. The resulting perturbative contribution obtained 
through Eqs. and (^, is given in Table || for different values of the 

strong coupling constant as(M^). In order to be conservative, and^ account 
for all possible sources of perturbative uncertainties, we have usedBLl 

= ±50A(‘‘)(as), (22) 

as an estimate of the theoretical error on Note that, for the relevant 

values of ««, this is of the saiiw size as K 3 thus, this error estimate is 

conservative enough to applyCl in the worst possible scenario, where the onset 
of the asymptotic behaviour of the perturbative series were already reached 
for n = 3,4. 

There have been attempts&il to improve the perturbative prediction by 
performing an all-order summation a certain class of higher-order corrections 
(the so-called ultraviolet renormalon chains). This can be accomplished us¬ 
ing exact large-n/ results and applying the naive non-abelianization prescrip- 
tion.o Unfortunately, the naive resuiJMiiajtion turns out to be renormalization- 
scheme dependent beyond one loop.E3’L3 More recently, a renormalization- 
scheme-invariant summation has been presented.eZI The final effect of the 








higher-order corrections (beyond K^) turns out to be small.| 


4 Power Corrections 


The 1/M^ contributions (5^^^ to the ratio Rt are simply the leading quark- 
mass corrections to the M^rbative QCD result of the previous section. These 
contributions are knownBcS to order a^. Quark-mass corrections are certainly 
tiny for the up and down quarks ~ —0.08%), byti^he correction from the 
strange-quark mass is important for strange decaysota {rriu = md = 0) 

^ulv/A = + ll.OSa^l, (23) 


where nris = ms{M‘^) is the running mass of the strange quark evaluated at 
the scale Mr- For as{M^) = 0.35, 6 uJ « —20%; nevertheless, because of the 
sin^ 9c suppression, the effect on the total ratio Rr is only —(1.0 ± 0.2)%. 

Since the r mass is a quite low energy scale, we should worry about possi¬ 
ble non-perturbative effects. In the framework of the OPE, the long-distance 
dynamics is absorbed into the vacuum matrix elements which are (at 

present) quantities to be fixed phenomenologically. If the logarithmic depen¬ 
dence of the Wilson coefficients on s is neglected (this is an effect 

of order a^), the contour integrals can be evaluated trivially using Cauchy’s 
residue theorem, and are non-zero only for I? = 2,4, 6 and 8. The corrections 
simplify even further if we also take the chiral limit (m^ = nid = rus = 0). The 
dimension-2 corrections then vanish because there are no operators of dimen¬ 
sion 2. In the chiral limit, sn(°)(s) = 0; thus only the 11*^°+^) (s) term in Eq. (||) 
contributes to Rr- The form of the kinematical factor multiplying 11^°+^^(s) 
in Eq. (^) is such that, when the s-dependence of the Wilson coefficients is 
ignored, only the D = 6 and D = 8 contributions survive the integration. The 
power corrections to Rt then reduce totra 


j( 6 ) 


j( 8 ) 

°ij,VIA 


~ -247r^ 




(D=6) 


M 6 






1 (£’= 8 ) 


Ml 


(24) 


and 0 for 2 n ^ 6,8. 

“ The present results do not include yet the resummation of the gn running corrections 
through the functions. 
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When the logarithmic dependence of the Wilson coefficients on s is taken 
into account, operators of dimensions other than 6 and 8 do contribute, but 
they are suppressed by two powers of as{M^). The largest power corrections 
to Rt then come from dimension -6 operators, which have no such suppression. 
Their size was first estimated in Ref. using published phenomenological fits 
to different sets of data: 


r(6) 


Re) 


(25) 


These power corrections are numerically very small, which is due to the fact 
that they fall off like the sixth power of Moreover, there is a large cancel¬ 

lation between the vector and axial-vector contributions to the total hadronic 
width (the operator with the largest Wilson coefficient contributes with oppo¬ 
site signs to the vector and axial-vector correlators, due to the 75 flip). Thus, 
the non-perturbative corrections to Rt are smaller than the corresponding con¬ 
tributions to Rt.v/a- a more detailed study of non-perturbative corrections, 
including the very small D = 4 contributions proportional to quark masses or 
to as(Af^)^, can be found in Ref. |[ 

The estimate ( p5| ) introduces a small uncertainty in the Rt predictions, 
since the actual evaluation of the non-perturbative contributions involves a 
mixture of experimental measurements and theoretical considerations, which 
are model-dependent to some extent. It is better to directly measure those 
contributpns from the r-decay data themselves. This information can be 
extractedE^ from the invariant-mass distribution of the final hadrons in r decay. 

Although the distributions themselves cannot be predicted at present, cer¬ 
tain weighted integrals of the hadronic spectral functions can be calculated in 
the same way as Rt- The analyticity properties of imply lEBS 



dsWis) 


1 

2 


/ 

J |s|=so 


dsW(s)n 


iJ) 

ij,V/A ’ 


(26) 


with W (s) an arbitrary weight function without singularities in the region 
s| < sq. Generally speaking, the accuracy of the theoretical predictions can 
be much worse than the one of Rt-, because non-perturbative effects are not 
necessarily suppressed. In fact, choosing an appropriate weight function, non- 
perturbative effects can even be made to dominate the final result. But this is 
precisely what makes these integrals interesting: they can be used to measure 
the parameters characterizing the non-perturbative dynamics. 

To perform an experimental analysis, it is convenient to use moments of 
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the directly measured invariant-mass 


distribution^ {k, I > 0) 



The factor (1 — s/sq)^ supplements (1 — s/M'^Y '^o Y in order to 
squeeze the integrand at the crossing of the positive real axis and, therefore, 
improves the reliability of the OPE analysis; moreover, for sq = it reduces 
the contribution from the tail of the distribution, which is badly defined exper¬ 
imentally. A combined fit of different R^iso) moments results in experimental 
values for as{M^) and for the coefficients of the inverse power corrections in 
the OPE. = Rt uses the overall normalization of the hadronic dis¬ 

tribution, while the ratios DY{M^) = R^(M^)/Rr are based on the shape of 
the s distribution and are more xlependent on non-perturbative efFects.E2l 

The predicted supfujessiontra of the^non-perturbative corrections has been 
confirmed by ALEPHB^ and CLEoJl3 using thejmoments (0,0), (1,0), (1,1), 
(1,2) and (1,3). The most recent ALEPH analysiscfl gives: 


(5np = ^ (cos^ + sin' dcSiY'^) = -(0-5 ± 0.7)%, (28) 

D>4 

in agreement with (|^). 


5 Phenomenology 

The QCD prediction for R^ is then completely dominated by the perturbative 
contribution non-perturbative effects being of the orderj^ the perturba¬ 
tive uncertainties from uncalculated higher-order corrections .LTH Furthermore, 
as shown in Table the result turns out to be very sensitive to the value of 
as{M^), allowing for an accurate determination of the fundamental QCD cou¬ 
pling. 

The experimental value for Rr can be obtained from the leptonic branching 
fractions or from the r lifetime. The average of those determinationsE2l 

= 3.649 ± 0.014, (29) 

corresponds to 

a^M^) = 0.355 ± 0.025 . (30) 

Once the running coupling constant as{s) is determined at the scale M^, it 
can be evolved to higher energies using the renormalization group. The size of 
its error bar scales roughly as a^, and it therefore shrinks as the scale increases. 
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Figure 2: Constraints on (5^^ and obtained from ALEPH dat 
the combined fit. All results are still preliminary. (Taken from Ref. 



The ellipse depicts 


Thus a modest precision in the determination of as at low energies results in a 
very high precision in the coupling constant at high energies. After evolution 
up to the scale Mz, the strong coupling constant in (^ 0 |) decreases toU 


a,(M|) = 0.1225 ±0.0025, 


(31) 


in excellent agreement with the direct measurementOtO at fj, = Mz, as{M^) = 
0.121 ± 0.003, and with a similar error bar. The comparison of these two 
determinations of in two extreme energy regimes, and Mz, provides a 
beautiful test of the predicted running of the QCD coupling. 

With as{M^) fixed to thej-mlue in Eq. (|^, the same theoretical frame¬ 
work gives definite predictionsu’LD for the semi-inclusive t decay widths 
Rt.a and Rt,s, in good agreement with the experimental measurements.aH 
The separate analysis of the vector and axial-vector contributions allows to 
investigate the associated non-perturbativemprrections {Rr.y — Rt,a is a pure 
non-perturbative quantity). Figure § showalj the (preliminary) constraints on 


Sy'^ and 5 ) 4 obtained from the most recent ALEPH analyses 
improvement over previous phenomenological determinations 

^ From a combined analysis of t data, ALEPH quotesE = 0.1225±0.0006exp ± 

0.0015th ± O.OOlOevol- 


?( 6 ) 


inary) cons 
lyse s RR .R A clear 
ionsoO is apparent. 
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Figure 3: Values of q:s(so) extracted^ from the R!^^{so) data. 0.0 The dashed line shows 
the three-loop QCD prediction for the running coupling constant. (Taken from Ref. 

The Cabibbo-suppressed width is very sensitive to the value of the 
strange quark massjj providing a direct and clean way of measuring iris- A very 
preliminary value, ms{M^) = (212 MeV, has been already presented at 

the last tau workshop.cS 

Using the measured invariant-mass distribution of the final hadrons, it 
is possible to evaluate the integral R^^{so), with an arbitrary upper limit of 
integration sq < M.^._The experimental sq dej^rdence agrees well with the 
theoretical predictions^ uputrpjather low valuesE3 of sq (> 0.7 GeV^). Equiv¬ 
alently, from the measuredBa distribution one obtains as(so) as a 

function of the scale sq- As shown£2l in Figure H, the result exhibits an im¬ 
pressive agreement with running predicted at three-loop order by QCD. 
It is important to realized that the theoretical prediction for does 

not contain inverse powers of Sq (as long as the s-dependence of the Wilson 
coefficients is ignored). The power corrections are suppressed by powers of 
1/M^; thus, they do not drive a break-down of the OPE. This could explain 
the surprisingly good agreement with the data for sq ~ 1 GeV^ 

A similar test was performed beforeEil for Rry, using the vector spectral 
function measured in e“'"e“ —> hadrons, and varying the value of the tau mass. 
This allows to study the behaviour of the OEE at lower scales. The theoretical 
predictions for Rry as function of agreeO well with the data for Mr > 1.2 
GeV. Below this value, higher-order inverse power corrections become very 
important and eventually generate the expected break-down of the expansion 
in powers of 1/M^. 
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6 Summary 


Because of its inclusive nature, the total hadronic width of the r can be rigor¬ 
ously computed within QCD. One only needs to study two-point correlation 
functions for the vector and axial-vector currents. As shown in Eq. (^, this 
information is only needed in the complex plane, away from the time-like axis; 
the dangerous region near the physical cut does not contribute at all to the 
result, because of the phase-space factor (1 — ■ The uncertainties of 

the theoretical predictions are then quite small. 

The ratio Rr is Jiery sensitive to the value of the strong coupling, and 
therefore can be useda to measure as{M^). This observation has triggered an 
ongoing effort to improve the knowledge of Rr from both the experimental 
and the theoretical sides. The fact that Mr is a quite low energy scale (i.e. 
that as{M^) is big), but still large enough to allow a perturbative analysis, 
makes Rr an ideal observable to determine the QCD coupling. Moreover, 
since the error of shrinks as /r increases, the good accuracy of the Rr 

determination of as{M^) implies a very precise value of as(M|). 

The theoretical analysis of Rr has reached a very mature level. Many 
different sources of possible perturbative and non-perturbative contributions 
have been analyzed in detail. A very detailed study of the associated uncer¬ 
tainties has been given in Ref. |^. The final theoretical uncertainty is small and 
has been adequately taken into account in the final as{M‘^) determination in 
Eq. (§). 

The comparison of the theoretical predictions with the experimental data 
shows a successful and consistent picture. The as{M^) determination is in ex¬ 
cellent agreement with the measurements at the Z-mass scale, providing clear 
evidence of the running of a^. Moreover, the analysis of the semi-inclusive 
components of the r hadronic width, Rr,v, Rr,A and and the invariant- 

mass distribution of the final decay products gives a further confirmation of the 
reliability of the theoretical framework, and allows to investigate other impor¬ 
tant QCD parameters such as the strange-quark mass or the non-perturbative 
vacuum condensates. 
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